ERROR ESTIMATES FOR BINOMIAL APPROXIMATIONS OF GAME PUT 

OPTIONS 



YONATAN IRON AND YURI KIFER 

INSTITUTE OF MATHEMATICS 
HEBREW UNIVERSITY 
JERUSALEM, ISRAEL 

Abstract. A game or Israeli option is an American style option where both the writer and the holder 
have the right to terminate the contract before the expiration time. As [9] shows the fair price for this 
option can be expressed as the value of a Dynkin game. In general, there are no explicit formulas for 
fair prices of American and game options and approximations are used for their computations. The 
paper 117) provides error estimates for binomial approximation of American put options and here we 
extend the approach of |17] in order to obtain error estimates for binomial approximations of game put 
options which is more complicated as it requires us to deal with two free boundaries corresponding to 
the writer and to the holder of the game option. 



1. Introduction 

A put option on a stock can be interpreted as a contract between a holder and a writer which allows 
the former to claim from the latter at an exercise time t the amount (K — S t ) + where K is a fixed amount 
called the option's strike, St is the stock price at time t and (x) + = max(x,0). In the American options 
case its holder has the right to choose any exercise time before the contract matures while in the game 
options case the contract writer also has the right to terminate it at any time before its maturity but 
then he is required to pay a cancellation fee in addition to the payoff above. 

The fair price of American options and of game options is defined as the minimal amount the writer 
needs to construct a self-financing portfolio which covers his obligation to pay according to the option's 
contract. It is well known that in the American options case the fair price can be obtained as a value of an 
appropriate optimal stopping problem while for game options we have to deal with an optimal stopping 
(Dynkin) game (see (9|). In general, both for American options and, even more so, for game options with 
finite maturity explicit formulas for their price are not available and approximation methods come into 
the picture while estimates of their errors become important. One of most easily implemented methods is 
the binomial approximation of stock prices modelled by the geometric Brownian motion and [17] provided 
corresponding error estimates for American put options. In the present paper we extend this approach 
in order to provide error estimates of binomial approximations for game put options. We observe that 
for perpetual game options some explicit formulas can be obtained (see [15 ]) but the finite maturity case 
studied here seems to be more realistic. 

Approximating the Brownian motion by appropriately normalized sums of Bernoulli random variables 
the paper |17] provided (error) estimates const-n -3 / 4 and const-n~ 2 / 3 for the difference between the 
price of an American put option and the price of its corresponding nth binomial model approximation. 
Using again the binomial approximation of the Brownian motion as above we construct in this paper two 
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approximating procedures such that the difference between the price of a game put option and its nth 
approximation in the first procedure is between const-n -3 / 4 and const-n -1 / 2 and in the second procedure 
is between const-n -1 / 2 and const-n -2 / 3 . The error estimates here are somewhat worse than in the case 
of American put options which is due to the lack of a smooth fit on the boundary of the writer's stopping 
region which causes substantial difficulties in the study of regularity of payoff functions. 

We observe that specific properties of game put options had to be used in order to obtain error estimates 
with the above precision. For instance, when payoffs are path dependent (and not only dependent on 
the present value of the stock) [lOj provides error estimates of similar binomial approximations only of 
order n -1 / 4 (m?i) 3 / 4 . Since price functions of game options can be represented as solutions of doubly 
reflected backward stochastic differential equations the results of [J] are also related to game options 
approximations. Nevertheless, approximations in [1] are not by binomial models, where computations 
can be done by means of the dynamical programming algorithm (see |10j), but by time discretizations, 
and so relevant probability space and c-algebras remain infinite which prevents effective computations. 
Furthermore, error estimates in [4] applied to our situation are of order n -1 / 4 , i.e. they are worse than 
for binomial approximations which we construct here for the specific case of game put options. 

Our exposition proceeds as follows. In Section [5] we provide basic results concerning game put option 
price functions, introduce our approximation processes and formulate our main result Theorem 12.11 In 
Section [3] we show that the price function can be represented as a solution of a variational inequality 
problem closely related to the Stefan problem (see [H]). We then use this representation to study 
regularity properties of the price function near the free boundary of the option's holder exercise region. 
In Section[3]we study the price function near the boundary of the exercise region of the writer. We use the 
information about this region from [14] in order to represent the price function as an explicit solution of 
the heat equation. This representation enables us to understand better the behavior of the price function 
near the boundary. We estimate also the rate of decay of the price function when the initial stock price 
tends to infinity. Section [5] is devoted to the proof of Theorem 12. II Finally, in Section [6] we exhibit some 
computations of the price functions and of the free boundaries. 

2. Preliminaries and main results 

The Black-Scholes (BS) model of a financial market consists of two assets among which one is nonrisky 
and the other one is risky. A nonrisky asset is called a bond and its price B t at time t is given by the 
formula B t — Boe rt where r is interpreted as the interest rate. A risky asset is called a stock and its 
price at time t is determined by a geometric Brownian motion 

(2.1) S t = S exp((r- y)* + «W t ) 

where n > is called volatility and Wt, t > is a standard Brownian motion defined on a complete 
probability space (fi, T, P). If So = % we write also Sf for St- The fair price of an American put option 
at time t with a strike (price) K and a maturity (horizon) time T < oo can now be written as a function 
F A {t, St) of time and the current stock price having the form (see, for instance, [13)). 

K 2 

(2.2) F A (t,x)= sup Eexp(-rr)(K -zexp((r- — )t + kW t )) 

reToT-t ^ 

where 7o, T-t denotes the set of all stopping times of the Brownian filtration with values in the interval 
[0,T— t] and E is the expectation with respect to the measure P. If we set ip (x) = (K — e x ) + , PA(t,x) = 
F A (t, e x ) and fi = r — ^- then we can rewrite (|2.2[) in the form 

(2.3) P A (t,x) = sup Eexp(-rr)V>(x + fir + nW T ). 

TdTo,T—t 

Relying on [9] (see also [15], [16] and [14]) we can also write the fair price of a game put option at 
time t with a strike price K, a maturity time T and a constant penalty 5 > as a function F(t, St) of 
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time and the current stock price in the form 

(2.4) F(t,x)= inf sup Eexp(-rcr A r)R(a, r) 

ae%,T-t T£T ,T-t 

where R(s,t) — (K — Sf) + + SI s<t and Iq is the indicator of an event Q. Using the functions P(t,x) = 
FA{t, e x ) and %j) as above we can rewrite this formula in the form 

(2.5) P(t,x)= inf sup E exp(— rcr A t) (ip(x + fur At + kW„ At ) + (5I cr<T ) . 

•?eTa.T-t TeTa, T -t 

It follows also (see [18], j9], [14], [16]) that the saddle point (optimal) stopping times for the game value 
expressions (I2.4|) and (|2.5p are given by 

(2.6) a* = inf {s < T — t : F(t + s,S*) = (K - S*) + + 5} A T and 

t* = inf{s <T-t: F(t + s, S*) = (K - Sf) + } A T. 
Next, we introduce our binomial approximations of the Brownian motion 

W t (n) = ^ V e k , t€[0,T], n = l,2,... 

where e^, fc = 1,2,... are independent indentically distributed (i.i.d.) random variables taking on values 
1 and -1 with probability 1/2 and [a] denotes the integral part of a number a. It is convenient to view 
{efcjfeLi as defined on the sequence space fi c = { — 1, 1} N = {£ = (£±,£2, ■•■) : £,i — ±1} by the formula 
e fc(£) = £fc if £ = (£i>£2,— )• Then W { will be defined on the probability space (S7 e , ^" £ , 7^ e ) where 
V t = {|, ^} N is the product measure and .F e is generated by cylinder sets. 

Now set 8* = Fa(0,K) which is the price of the American put option with a maturity T and a 
strike A. Clearly, if the penalty 8 > 8* then it does not make sense for the writer to cancel the 
corresponding game put option, and so in this case the prices of American and game options are the 
same, i.e. Fa{0,K) = A(0,A). Since approximations of American options were studied in |17] we 
assume in this paper that 8 < 8*. Observe that FA{t,K) is continuous in t and it is strictly decreasing 
to as t increases to T, and so for each 8 £ [0, 8*} there exists a unique tg < T such that FA(ts, A) = 8. 
Furthermore, we can define k n to be the minimal k <E N such that 8 > FA(Tk/n,K) and set /3^™- > = Ssa, 
In order to define two sequences of functions Pj and P% , n = 1, 2, ... which will approximate P(0, x) 
we set Xj: — x + and introduce stopping times 

(2.7) ct (,i) =inf{t € [0,f3 {n) ) : In A - \n\h - 2nVh < fJ,h[U + x[ n) < In A' + \fi\h + 2kVJi} A T 

where a^ n ' — T if the infimum above is taken over the empty set. Introduce a filtration {Q t = F[ t /h] , t > 0} 
where is the trivial c-algebra and JFk is generated by ei, e k - Denote by the set of all stopping 
times with respect to the filtration {Gt} taking on value in the set {kh, k = 0, 1, ...,n}. Then, clearly, 
a (n) g q-(n)_ NoWj for x <\ n K we define 

(2.8) P[ n \x) = sup E(e-™ ( " ,AT (^T + 4»))l {r < ffW}+ fl {ffW<T} )) 

r£T<") 

and for x > In A we set 

(2.9) P 1 (n) (^)= sup E(e-- < '* ) ^(VXMT + A("))I {T < ff( „ )} + (5-A e ^l^+^)l {CT( „ )< ^ 

rer<") 

The second approximation function is defined for all x by 

(2.10) P 2 (n) (x)= sup Ete^'^/xT + A^^ 

We can formulate now our main result. 
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2.1. Theorem. For each x there exists C = C{x) such that 

(2.11) - JL< P («) (;r) _p ( o, x )<-g_ and -^<pW (l )-P(0,,)< 



In the following sections we will analyze regularity properties of the price function P(t, x) of game put 
options and will complete the proof of Theorem 12. II in Section [S] providing some computations in Section 
[HJ The general strategy of the proof resembles that of [T7] but the study of the price function of game 
put options is more complicated than in the American options case, in particular, because of appearance 
of two exercise boundaries (holder's and writer's) having different properties. Our proof will be based on 
regularity properties of solutions of parabolic partial diferential equations with free boundary and of the 
corresponding variational inequalities and we will rely also on some prior results from [17] , |15) and |14j . 

3. Price function near the holder's exercise boundary 

3.1. Some previous results. First, we state the following result from [14] (see also [16]) which we will 
use later on. 

3.1. Proposition, (i) There exists an increasing function b : [0,T) — ► R such that lim^T b(t) = K 
and F(t, x) = K — x for all (t, x) satisfying < x < b(t). 

(ii) There exists < 5* such that for every < 5 < S* there is a — (5(5) > so that F(t, K) = 8 for 
t € [0, 0] and for t > we have F(t, x) = F^(t, x) for all x > 0. 
(Hi) Furthermore, 

k 2 

F t {t, x) + x 2 F xx (t, x) + (r - —)xF x {t, x) - rF(t, x) = 

for all 

(t, x) € (0, T) x R+ \ ({(*, x) : x < b(t)} U [0, 13) x {K}). 

In particular, F(t,x) is of class C 1,2 , i.e. continuously differ entiable once with respect to t and twice 
with respect to x, and so, in fact, it is a smooth function there. 

(iv) Finally, F(t,x) is convex and strictly decreasing in x and nonincreasing in t. 

Next, we introduce an operator T> which acts on Borel functions u(t, x) on [0, T] x R by 

(3.1) T>u(t, x) — — [u(t + h,x + K\/h) + u(t + h,x — nVh)] — u(t, x) 

Clearly, \V{t,x) can be viewed as a discretization of the differential operator + ■ We will rely 

on the following results from |17j concerning the operator T>. 

3.2. Proposition. For each Borel function u on [0,T] x R there exists a martingale (M t )o<t<T with 
respect to the filtration Q t , t > such that Mq = and for every t £ {0, h, 2h, T}, 

nt/T 

(3.2) u(t, X[ n) ) = u(0, x) + M t + Vu((j - l)h, X$l 1)h ). 

3 = 1 

3.3. Proposition. Let < t < T — h and x £ R. Assume that u is a C 2 function on ([t,t + h] x [x + 
Kyh, x — K-yh] ) . Then 

(3.3) Vu(t,x) = - dy dz(z-—{t + y 2 ,x + z) + 5(u)(t + y 2 1 x + z)) 

K Jo J-kv otdx 

where 

k 2 

5(u)(t, x) = u t (t, x) + —u xx (t, x). 
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We will need also the following result concerning the free boundary s(t) = \n(b(t)) of the holder exercise 
region of our game put option which in the case of American options appears as Proposition 1 in |17j 
and it can be proved for game options in the same way. 

3.4. Proposition. Let < ti < t 2 < and let x = s(0) < x\ = s((3) < \nK then (s(ti) - s(t 2 )) 2 < 

\P(t 1 ,x)-P(t 2 ,x)\. 

We also observe that it follows from the Berry-Esseen estimate that for some constant C\ > inde- 
pendent of j, n > 1, 

(3.4) P{\X$-z\< K Vh}<^. 

We will also rely on the following standard bounds on derivatives of solutions of 2nd order parabolic 
equations with constant coefficients (see, for instance, [3] and [5]). 

3.5. Proposition. Let D — (0,T) x (0,1) and let w(t,x) £ C[D] be a solution in D of the following 
parabolic equation 

k 2 

—w xx + fiw x -rw = w t . 

Suppose that w(0,x) = for all < x < 1 and that there exists A > such that \w(t, x)\ < A for all 
(x, t) £ D. Then for every fc, n and < a < b < 1 there exists C = C(k, n, a, b, T, A) such that 

d k+n w 

(3.5) \——( t ,x)\<Cforall(t,x) £ (0,T) x [a, 6]. 

3.2. Price function and variational inequalities. Next, we will show that the price function of the 
game put option can be represented as a solution of a variational inequality (v.i.) problem which is a 
generalization of the Stefan problem (see pj] ,VIII). This will enable us to derive certain regularity 
properties of this price function which we will use later on. Details of some of the proofs concerning the 
solutions of the v.i. problem below which are similar to the proofs in the case of the Stefan problem will 
not be given here. For the corresponding results in the American put option case we refer the reader to 
[13], [17] and to references there. 

Let T be such that j3 < T < T and set 

k 2 d 2 d k 2 

(3.6) A = -— + „_ - r where M = r - -. 

Using the maximum principle, properties of price functions of American and game put options and the 
fact that after time f3 the price functions of the game and American option are the same we obtain that 
for every x > s(T') the time derivative Pt(T',x) = PA,t{T',x) is strictly negative and we can find a, b 
satisfying s(T') < a < b < In AT such that for some constant c > 0, 

(3.7) -P t (T',x)>c VxG[a,6]. 

Relying on Proposition 13. If hi) we also observe that for all (t,x) £ [0,T"] x (s(t), In K), 

(3.8) ^(t,x)+AP(t,x) =0, P(t,x)>K-e x V(t,x) £ [0,T'] X (s(t),1nK), 

P{t, x)=K~e x yt£ [0, T'}, Vx < s(t) and P t < 0. 

Let ao be such that an. < s(0) < s(T r ) < b. Introduce the domain D = (0, T") x (ao, b) and for all (t, x) 
in the closure D of D define the functions 

(3.9) v(t, x) = P(T' -t,x)- P(T', x) and f{x) = AP(T', x). 
We obtain that 

(3 10) f(x)-i - p t( T '> x ^ s(T')<x<b 

[6Ai)) TW-i _ rK> a <x<s(T') 



6 



Y.Iron and Y.Kifcr 



and from the definition of v(t, x) it follows that for any (t, x) € D, 

(3.11) P t (t, x) = -v t (t, x), P tx {t,x) = -v tx (t,x), P tt (t,x) = -v tt (t, x), 
P x (t,x) +P x (T',x) = -v x (t,x) and P xx (t,x) + P xx (T',x) = -v xx (t,x). 

Since P x (T',x) and P xx (T',x) are bounded we obtain that the integrability properties of the first and 
second order derivatives of P(t,x) and v(t,x) are the same in D. Now set 

(3.12) ip(t)=v(t,b) , g(t) = v t (t,b) for < t < T'. 
Then by (02]) and (j3~TT|) . 



(3.i3) m 



f g(r)dr = f v t (r, b)dr > for < t < T' . 
Jo Jo 



It follows from $3Jj$ and ([3~9l) - ([3~T0| that on the set v > 0, 

(3.14) v t - Av- f = -P t (T' - t, x) - AP(T' -t,x) + AP(T', x) - f(x) = 
and on the set v — we obtain 

(3.15) v t - Av - f = rK > 0. 
Hence we arrive at the following (see [TTj'). 

3.6. Lemma. The function v is the unique solution of the following variational inequality problem, 
v.i. Problem 1: Find v G L 2 [0, T'\ H 2 (a , b)] n H X [D] such that 

(i) v,v t > 0. 

(ii) (v t — Av)(w — v) > f(w — v) a.s for every w € L 2 [D], w > 0. 
(Hi) v(t, b) =%/) for <t <T', x = b. 

(iv) v(t, ao) = for < t < T' , x — oq. 

(v) v(0, x) = for t = 0, a <x<b. 

Proof. We shall prove uniqueness, the fact that v is a solution to v.i. Problem 1 follows from (|3.9|) - (|3.15|) . 
Assume that v and v are two solutions of v.i. Problem 1. Since v > (property (i)) we can use the 
property (ii) of v and replace w by v. Since both of them are solutions we obtain that 

(3.16) (vt — Av)(v — v) > f(v — v) and (vt — Av)(v — v) > f(v — v). 

Define the parabolic boundary as the boundary of D without the interval {T'} x (ao, b) and let u = v — v. 
Note that u is zero on the parabolic boundary and the sum of the two inequalities (|3.16p is 

(3.17) u t u — —u xx u — fxu x ii + ru 2 — (u t — Au)u < 0. 

Integrating both sides of p,17p on (0,T') x (ao,b) we obtain four terms on the left side. For the first 
term we have 

i-b rT' pb ^ 

u(t,x)u t (t,x)dtdx = / — u 2 (T', x)dx > 0. 

ao JO J ao ^ 

Integration by parts of the second term and the fact that u = on the parabolic boundary yields 

K 2 f T ' f b K 2 f T ' f b 

(3.18) — / / u xx {t 1 x)u{t,x)dxdt = — I I u x {t 1 x)dxdt > 0. 

* JO Ja ^ JO Ja 

For the third term note that u x u — bt 4 - and that u(t, ao) = u(t, bo) = for every t, and so 

I u x (t,x)u(t,x)dxdt = 0. 

JO Ja 

The last term satisfies r J a<j u 2 (t, x)dxdt > since r > 0. We conclude that the left side of f|3 . 1 T[) can 
not be negative and so it must be zero. Since all terms in the left hand side of (|3.17p are non-negative 
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and their sum is equal to we obtain that r J aa u 2 (t, x)dxdt = 0, and so u = almost everywhere 
(a.e.). Hence, v = v a.e., and so there is only one continuous solution. □ 

Denote parts of the boundary of D = (0, T") x (do, b) by 

T 1 = [0,T r \x{b}, r 2 ={0} x (oo,6), r 3 = [0,T'] x {a }, T = Ti U T 2 U T 3 

and set 

9 d k 2 d 2 d 

Thus, r is a parabolic boundary of D. For every e > we define following functions. 

(1) A smooth function / (e) (a;) > f(x) on (a ,6) such that / (e) (a;) = /(a;) for s(T') < a < x < b and 
for a < x < ai where a\ satisfies ao < a\ < s(T') and lim e ^ f^ e '( x ) = f° r a o — x — 

(2) A smooth function /3^ e \v) satisfying 

/3 (£ \v) = for all v > e, /3 (e) (0) = -1, ^{v) > and (v) < 0. 

(3) ^ e \t) = ijj(t) + e with V defined in ([3l2]l . 

(4) A smooth function 77(3;) such that < 77(2;) < 1 and for some a < a 2 < b, 

f](x) = 1 for d2 < x <b and 77(21) = for ao < x < a. 

Set F e {x,v) = f ie) (x) - rK0&(v) which is a Lipschitz continuous function and for every constant C 
there is Mq such that C\F^ (x,v)\ < M whenever M > Mq and |v| < M. Let <j)^ be a function on V 
satisfying 

^\ Tl = ^)(t), ^\ V2 = er,(x), ^\r 3 = 0, 
and, moreover, relying on Chapter 3 in [5] we can choose 4>^ so that 

(1) <f>^ € Ci+s\D] for some < 8 < 1 (in fact for each 5) and we refer the reader to Chapter 3 in 
[5] for the definition of C2+s[D] and for conditions yielding that a function defined only on the 
boundary T can be extended to a function from C2+5 [D] . 

(2) = F^ e \x,ip) at the points (0,6) and (0,a ). 

By the theory of semi- linear parabolic equations (see [5]) there exist a function € C^+^-D] for some 
< 7 < 1 such that 

(3.20) Ltj( £ ) = (x, v {£) ) and v (e) | r = </> (e) . 

In particular , , i4x , , are continuous on 13. 

Let to = Tjj 6 * 1 . By differentiating with respect to t the equation (|3.20[) and taking into account f|3.12[) , 
(13.201) and the properties of cj)^ we obtain that 

(3.21) w t - ^w xx - nw x + (r + rK$\vM))w = 

where w(t, b) = g(t) V0 < t < V and w(t, o ) = V0 < f < T" 
w(0,x) = /^(x)+e(4^(a;)+Wx(x)-r77(x)) -r^M e) (7j( £ )(0,x)) Va < z < 6. 

We see that in D the function id is a solution to a parabolic equation and since r + > we can 

use the maximum principle 

(3.22) min(min(u> ),0) < w(t,x) < max(max(7i;),0) V(*,x) e D. 

Therefore in order to bound the function w we only need to bound its values on the parabolic boundary. 
First, we estimate the left hand side of p.22j) . For a < x < b we have that v^(0,x) = srj(x) < e, and so 
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0& ( V (s)) < 0. In view of p.7p . (|3.10l) and the definition of f( e ) above there exists £q > such that for 
every < e < eq, 

w(0, x) > f(x) + e(r) xx (x) + /j,r) x (x) - rr)(x)) 
= -P f (T' ,x) + e(^r) xx (x) + nr/ x (x) - rr](x)) >0 V a < x < b. 

On the interval ao < x < a we have r\ = and since v^ e )(0,x) = er](x) =0we see that (i/ £ ) (0, x)) — 
— 1. Since on this interval f{x) > —rK we obtain 

w(0,x) = / (e) (x) - rK0 {£) (0) > f(x) + rK > for a < x < a. 

Hence, w > on T2- We obtain next that, 

dP 

w{t, b) = —q^{T ~ t, b) > on Ti and w(t, b) = on T 3 . 

It follows that min (minr(w), 0) = 0. 

Next, we estimate the right hand side of (I3.22[) . On I^ we have that 

w(0,x) < \f (£ Hx)+e(^-7 lxx (x)+ mx (x)~rr ] (x))~rK0^(v^(O,x))\ 



2 
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< sup \f (£> | + e sup | —r) xx + (j,r] x - rrj\ + rK < C 
where C > is a constant independent of e, and so 

= min ( min(w ),0) <w{t,x) <max(max(- — (T' - t,b),C )) =C X . 
We conclude that there are some constants Eq and C\ such that for every < e < £q, 

(3.23) < v ( t E \t,x) < C x . 

Since v[ e ^ > and v^ e ) (0, x) > for ao < x < b we deduce that i/ 6 - 1 (t, x) > and because uf' is uniformly 
bounded it follows that y( e ) is also uniformly bounded. By the properties of 0( e ) we see that 

(3.24) - 1 < (5 {€ \v (e) ) < or \\0 (e) {v (e) )\\ L ~ [D] < 1. 

Let Dq — (0, T) x (d2, &) be an upper subrectangle of D where a2 is the same as in the definition of the 
function 77 in (4). From the definition we have i)( E )(0,x) = Er]{x) = e in Dq and since u E) j is nonnegative, 
we obtain that v^(t, x) > e, and so {e) (w (e) (t, x)) = 0. 

This means that on Do the function v^) satisfies the parabolic equation 

For w = v\ £ ^ and < e < £q we also have that 

Lw(t,x) = V(f,x) € L>o and uj(0,x) = / (e) (x) when a 2 < x < b. 

Next, let y(t,x) be a function on Do such that 

Ly(t,x) = W(t,x) e D , y(0,x) = f (e) {x) Va 2 <x<b 

and all of its first and second order derivatives are bounded there. Such a function exists since we can 
choose a smooth function on the remaining part T \ {0} x [a 2 ,b] of the parabolic boundary r of D 
which extends f( e >(x) as a smooth function to the whole Do, and then use Theorem 12 from Chapter 3 
in For each e < £0 we define z[t, x) — w(t, x) — y(t, x) in the domain Do where w(t, x) = v\ (t, x). 
Then z(0, x) — w{Q,x) — y(Q,x) — for every a 2 < x < b. Fix xo G (&2,6), then by Proposition 4.5 
from Section 4.1 of [IT] we obtain that \zt(t,Xo)\, \ztt(t,xo)\ < C for every < t < T' where a constant 
C > is independent of e. Since we assume that \yt(t, xq)\, \ytt(t, xo)\ < C\ for < t < T' it follows that 
\w t {t, Xq)\, \w t t(t, xq)\ < C\ + C (for every e). Let D\ = (x ,b) x (0,T), then by Theorem 6 in chapter 3 



Approximations of game put options 



of [5] we obtain that for every e > (and in fact every < a < 1) v^ 6 \ v\ e ' G C2+q:[-Di] and there is 
constant C independent of e such that 



\^ ) k +a + K > \ 2+a <c. 



In particular, we get 

(3-25) H e) \\ L ^ iDl) + \\vg\\ L ~ w + ||i£ e) |U~ (Dl) < C. 

Considering again the whole region D we have 



2 * 



(t, y)d tf = - / (M £) (t, 2/) " v { t £) (t, y) - rv^ (t, y) + ^ (v^ (t, y)) + (y)) dy. 



Hence, 

v x s) (t, x) = vi e) (t, b) + £ [»(v^ (t, b) - v& (t, x)) 
+ J b (- v[ e) (t, y) - rvOO {t , y) + 0& («<«> (t, y)) + f& (y)) dy] . 

Since all terms in the right hand side are uniformly bounded there is a constant C > independent of e 
such that ||i>x || < C for every < e < £o- Now we see that in the equation 

4} (*, y) = 4 (*, i/) - « t (£) (*, v) - ™ (£) (*, y) + P {£) (*, v)) + f (e) (y) 

all terms in the right hand side are uniformly bounded and therefore the term in the left is uniformly 
bounded, as well. 

We summarize this in the following lemma. 

3.7. Lemma. There are constants C > 0,£o > such that for every e < Eq, 

Hlh-ID] + H e) \\L- lD] + ||w t (B) |U-[l,] + \\V^\\ L « [D] < C. 

We now obtain the following (see [TTj). 

3.8. Proposition. For any 1 < p < oo and t g [0, T'], — > v as e — > weakly in W 1,P {D). Further- 
more, — > v uniformly on D and also v x — > v x uniformly in x € [0, if] for each t € [0, T']. The 
function v is the unique solution of v.i. Problem 1. 

Next, we analyze properties of second order derivatives starting with the following result. 

3.9. Lemma. There is a constant C > such that for any < e < Eq, 

(v { t e J(t,x)) 2 dxdt < C. 

J a 

Proof. Set v = v (<0, /3 = /3( £ ) and w = v t e . Multiply the equation (|3.21l) by w to obtain 

wwt —ww xx — fiww x + (r + rKfi 1 (v))w 2 = 0. 

Integrating this equation over (ao,b) and recalling that j3'(v),t and K are non- negative we obtain that 
for any < t < T', 

1 d f b k 2 f b f b 1 dw 2 

(3.26) -— / w 2 (t,x)dx — / w(t,x)w xx (t,x)dx — fi / -— — (t, x)dx < 0. 

2 at ./„„ 2 ./„„ ./„„ 2 da; 



By (pHUl) and (EH3")) we estimate the third term in (pHBf 
<- b 1 dw 2 



A* / Y^-{t,x)da\ = \v\\w 2 {t,b)-w 2 (t,a )\ = \^\w 2 (t,b) < C 2 . 
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For the second term in (I3.26|) we see that 

— o~ / w(t,x)w xx (t,x)dx = — ( / wl(t,x)dx-w(t,b)w x (t,b)+w(t,ao)w x (t,a )). 



Since w(t, ao) — and the function w(t, x) is uniformly bounded in D we see that w x — v\ E ] is uniformly 
bounded near the boundary [0,T'] x {b} and w(t,ao)w x (t 7 a ) = while \w(t,b)w x (t,b)\ < Ci for some 
constant C2 > independent of e. Thus, we conclude from (|3.26l) that 

1 d f b k 2 f b 

Ydt] w2 ( t ^ x ) dx + Yj w l(^ x ) dx - C 3 

for some C3 > independent of e. Integrating the last equation over [0, T'\ we obtain 

k 2 r T ' r b 1 r b 

— I / w 2 x {t,x)dxdt + - I (w 2 {T,x) - w 2 (0,x))dx < C 3 . 

2 Jo Ja 2 J ao 

Since the function w is uniformly bounded it follows that there is C > independent of e such that 

rT' rb 

w 2 (t, x)dxdt < C. 

J ao 

□ 



We will now deal with the L 2 properties of the function v^\t, x). 
3.10. Lemma. There is a constant C > such that for any < e < Eq and every < a < t < T' . 

"\v£\t,x)) 2 dx+ f I {v\f{x,t)fdxdr< C 



Proof. Set v = v (0,/3 = pi*) and w = vf' . Multiplying (|3.21[) by the function wt we have 

k 2 

wf —w xx w t - fiw x w t + (r + rK/3'(v))ww t = 

and an integration with respect to x over (ao, 6) yields 

f b k 2 f b f b f b 

(3.27) / w 2 dx — / WxxWtdx — fi / w x Wfdx + / (r + rK(3'(v))wwtdx = 0. 



Fix some t € [0, T]. Since uit(t, ao) = w(t, ao) = we see that 

k 2 f b k 2 k 2 d f b 

w xx {t,x)w t (t,x)dx = —w x (t,b)w t (t,b) - —— J w x {t,x) 2 dx. 

From (|3.25p it follows that ^-\w x (t,b)wt{t,b)\ < C\ for some constant C\ > independent of e, and so 
we obtain 

(3.28) J b o w 2 x (t, x)dx + J b o w 2 dx + J b Jr + rKp'{y))w(t, x)w t {t, x)dx 

<H J b g w x (t,x)w t (t, x)dx + Ci. 
Now we deal with the last term in Q3.27p . Since /3"(v) < and v,w>0we obtain that 
fjr + rKp{v))ww t dx - \ f b Jr + rKp'{v))f t w 2 dx 
= hi fa (r + rKf3'(v))w 2 (t,x)dx - \ £ rKp»{v)wHx > \fj b jr + rK (3' (v))w 2 dx. 
We plug this inequality into Q3.28[) and obtain 

Id f b k 2 f b f b 

— — / [— w x (t, x) + (r + rK (3 1 (v))w 2 (t, x)]dx + / w 2 (t,x)dx+ < /i / w x (t, x)wt(t, x)dx + C\. 
2 dt J ag 2 J ao J ag 
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Integrate the last inequality with respect to r' over the interval (t, t) to obtain 

I C *) + + rK0>(v))w*(t, x)]dx + J* /* u£(t, *)aW 

< m£ lac M*>*)IM*. + d(i - r) + § /^[^(r, x) + (r + rf/3'(«)) w 2 (r, i)]dx. 

Next, integrating in t over the interval (0,cx) for some < a < t and taking into account that (r + 
rKf3'(v))w 2 > by the property (2) of /3 we obtain that 

(3.29) § j£ ^(i.^ + ^^/^Ki.rrMxdr'dr 

<C 2 + \n\S°£? aQ K(t,a:)|K(t,x)|dxdT'dr + ^^/^[^(r.a;) + (r + rK(3'(v))w 2 (r, x)]dxdr. 

Now, by (|3.23p , (|3.24p and Lemma 13.91 together with the Cauchy-Schwarz inequality we estimate the 
right hand side of (|3.29[) by a constant C3 > independent of e. Hence, 

CK 2 f a 0~K, 2 

C3 > — — / w 2 dx + / / w 2 dxdr' dr > — — / w 2 dx + a I I w 2 dxdr' 

^ J do « J t Jflfl J ao J a J aa 

and Lemma [3.101 follows . □ 
As a corollary of previous results we obtain 

3.11. Proposition. Let (3 < a < T and a < s(0) < s(cr) < b < InK. Define D a = (0, a) x (a, &). T/ien 

(3.30) P{t,x) e H 2 (D a ) 

where by definition H 2 [U] is the set of all the functions in L 2 [U] with an L 2 weak second order derivatives. 
Also there exists C > such that for every < t < T' , 

(3.31) /jS(^)l^=llS(^)II^M]<^ 

Proof. From Lemma [3.101 Lemma 13.91 and Lemma 13. 71 we obtain that {v^} e<eo are uniformly bounded 
in H 2 [D a ] and so they have a weak limit v € i7 2 [Z? cr ]. Since u' e ) — >• v uniformly we must have that v = v, 
and so v € H 2 [D a ]. Since v is the solution of (|3.15|) we can apply Proposition 13.51 and using the fact that 
the constant C in (|3.5|) doesn't depend on t we can obtain in a similar way that for a fixed ct there is a 
constant C > such that for every < t < a, 

IK,t(*, -)IU 2 [ao,fc] < C- 

From p. lip we can deduce the same result for the function P(t,x). □ 

3.12. Corollary. For each < t < T the function v t {t,x) is Holder continuous with a Holder exponent 
1 

2 • 

Proof. For every < t < T Proposition 13.111 gives us that v t (t,x) S if 1 [an, 6]. Hence, the result is a 
consequence of the Sobolev inequality. □ 

3.13. Corollary. For every < t < T' the functions Pt(t,x) and P xx (t,x) as functions of x are continuous 
in the closed interval [s (£),&]. 

Proof. For the function P t (t,x) the result follows from p. lip and the previous corollary. Since P(t,x) is 
a solution of (|3.8p in the interval {(t,x) : s(t) < x < In if} and since the functions P x (t,x) and P(t,x) 
are continuous in the interval [s(t), b] we obtain the result for P xx , as well. □ 

3.14. Corollary. Let [3 < a < T and a < s(0) < s{a) < b < InK. Define E a = {(t,x) : < t < 
f3, a — fit < x < b — fit} and u(t, x) — e~ rt P(t, x + fit). Then 

(3.32) u(t, x) G L 2 [E a ] 



12 



Y.Iron and Y.Kifcr 



and there exists C > such that for every < t < ft, 

(3.33) J^\ — { t,x)? d x<C. 

Proof. The assertion (13.321) follows from Proposition 13.111 and the definition of u(t,x). For (|3.33|) note 
that 

d 2 u , , rt/ dP . . a 2 p, , a 2 p, NN 

^ (*,*) = e"* ( - r- (t,x + 0) + (t,x + Mt) + ^ (t, * + *)) , 

then use (|3.31|) and the fact that for (t,x) G E a the functions j^(t,x + /it) and j^(t,x + fit) are 
bounded. □ 

4. Price function near the writer's exercise boundary 

4.1. Regularity properties of price function. Let F(t,x) be the price function of the put game 
option (see Section [2]). We begin this section by showing that near the writer's exercise region Fi = 
{(t, K) : < t < ft} the function ^ is continuous. Let 

(4.1) Y} s - X] = (Y t Us ' x \ Y t 2 ' M ) = (s + t, Sf) 

which is a non homogeneous in time Markov process in R + x R where Sf = X e^ t+KBt and \i = r — 
Let 

, , . d k 2 x 2 d 2 d 

(42) LY =dt + —dx^ + rX 8x- r 

which is the infinitesimal generator of Y t when considered on the space of all C 2 functions. This is a 
parabolic operator with bounded smooth coefficients in the domain 

(4.3) D=(0,ft)x(k,K) 

where k > 0. Let Pr S)X ] and E[ s ,x] be the probability and the corresponding expectation for the Markov 
process Y starting at the point [s, x\. We will first show that for any to G [0, ft), 

(4.4) lim P [t , x] [y T er 1 ] = i 

(t,x)->(t ,K) 

where t = t(T) and for any closed set Q C R+ x R we set t(Q) to be the arrival time at the set Q 
for a Markov process under consideration which is Y t here. Indeed, choosing an appropriate nonnegative 
function <f> < 1 on the boundary T and relying on Chapter 3 in [5] we can choose u(t, x) G C 1,2 (D) which 
solves the equation Lyii = in D and equals 1 on the boundary part Ti for < t < t\ < ft while decaying 
smoothly to when t grows to ft. Then 

u(t,x) = E M 0(r T ) < p m {f t g ri}, 

and so 

1> liminf Pr t X ]\Y T G T{\ > lim u(t, x) = u(t , K) = 1. 

(M)-Kto,h) ' (t,x)->(to,K) 

Next let f(x) = (K — x) + and g(x) = f(x) + 5. Recall that the price of a put game option with an 
expiration time T and a constant penalty 5 can be written in the form 

F(t, x) = sup inf _ Ju x] (/, g, a, r) 

(1<t<T»<"<T 

where T = inf{t : Y t x — T} and for any bounded Borel functions / and g we write 
J [m1 (/,3,<t,t) = E [M ,[e-™ AT ( S (}; 2 )l {ff<l} + /(y T 2 )I {l < 5} )]. 

Set 

(4.5) f s {x) = F(s,x) when ft < s < T and F s (t,x) = sup a < T < s info< a < s J [t ^ x] (f s , g,a,r). 
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where s = inf{w : yj 1 ^ — s}. Let < cr*,r* > and < cr*,r* > be the two saddle points (see [5]) 
corresponding to the optimal stopping games with values F(t,x) and F s (t,x), respectively, and so 

(4.6) a* = inf{0 < t < f : F(Y t ) = g(Y t 2 )} r* = inf{0 < t < f : F(Y t ) = f(Y t 2 )} 
a* 8 = inf{0 < t < S : F s (Y t ) = g(Y 2 )} t* = inf{0 < t < S : F s (Y t ) = f s {Y t 2 )}. 

Then 

F(t,x) = J [t>x ] (f,g, <t*,t*) and F s (t,x) = J^ x] (f s ,g,cr*,T*). 

4.1. Lemma. For all < t < s < T and x > 0, F s (t,x) = F(t,x). 
Proof. We have 

F s (t,x) = J [ttX ](f s ,g,a*,T*) < J [t , x] (f s ,g,a*,T*) 

= E [t)X] [e- rff * A < (f s (Y 2 )l {T ,< a * } + F(y cr .)V<r;})] 
< E M [e- r » ,AT :(F(y T; )i {T; < ff . } + i?(y tr .)V<r;})] = v [t ,4e~™ * A <F(y ff . AT .)] < F(t,aO. 

Indeed, the first inequality above follows by the saddle point property. The second inequality holds 
true since F is nonincreasing in the time variable, t* < s = s — t for Y^' x ' and Y 1 '^'^^*) < s. The 
third inequality is satisfied since the process e~ rY °* A ™F(Y a * Au ) is a continuous supermartingale in u with 
respect to ~P\t, x ] ( see [S])- For the other direction we have 

F(t,x) < V [t , x] [e- rSA °>T 'f (y- AffJAT .)] = E [t)X] [ e -" ACT .* A -* (/(y T 2 , )I T -<^ 

+F(y s )I s<T . ACT , + 5 (y £T 2 ; )I CT;<SAT -)] < E^Ie-'-^^^/.O^^JIr-A^. 

+ 5 (y ff 2 ; )I CT . <SAT -)] = J [t)X] (/ s ,<?,<,T* AS)) < J [M] (/ S , 5 ,a s *,r;) =F s (t,x) 

where we use the submartingale property of e~ rY ^*^F{Y T */\ u ) in u. □ 
Now for any bounded Borcl functions / and g set 

I s (t,x,f,g)= sup m£ J[ tx ](f,g,cr,T). 

0<T<S°< a < S 

From the time homogeneity of the process Y 2 = St we obtain that 

(4.7) I s+h (t + h 7 xJ,g)=I s {t,x,f,g). 

4.2. Proposition. There is a constant C > such that for any (t,x) £ (0,(3) x (k,K), 

dF 

o<- ¥ (t,x)<cp M [r;< ff : +)l ]. 

Proof. The left hand side of the above inequality follows from (iii) and (iv) of Proposition 13.11 For the 
right hand side, let h > Q be such that (3 + h < T - h and t + h < (3 and let j3 < s < T — h. By (see [L]) 
the price function of an American put option has a bounded derivative with respect to t in [0, s + h] xl, 
i.e. C — sup( tiX ) g r 0)S+fc i ><R 1 9Fj dt ^ I < 00 ' T ms together with Proposition 13 . 1 f ii) yields 

(4.8) sup |^M|<C. 
Next, by Lemma |4. H and the saddle point property, 

(4-9) F(t,x) = F s (t,x) = J [t ^(fs,g,<J* s ,T* s ) < J [ttX] (f s ,g,a* s+h ,T*). 

By Lemma \4. 11 (14.71) and the saddle point property, 

(4.10) F(t + h,x) = F s+h (t + h,x) = I s+h (t + h,x,f s+h ,g) = I s (t,x, f s+h , g) 

= J[t,x]{fs+h, 9,0-* s+h ,T* +h ) < J[t, x ](f s+ h,9,0-* +h ,T*). 
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Now, gH) , (gU) and (I4TTU1) yields that 

< ^[e^X^^^j] < CP [M] [r* < 
Passing to the limit as h — »• we obtain the result. □ 
4.3. Corollary. For every < to < (3, lim^^^K) %{t,x) = 0, and so lim^ ^j — ^(to,ln k) ^f^^) = 0. 
Proof. In view of Proposition 14.21 we only have to show that for every < to < f3, 

(t,x)-+(t ,K) 

Let D be as in (|£3| . T 2 = {(P,x) : k < x < K} and T 3 = {(t,k) : < < < 0}. It follows from the 
definition of t* and o-* s+h that for every (x, i) G £>, 

{r(ri) < r(r 2 n r 3 )} C {o-* +h < t*} with respect to P[*, x ]- 

From (|4.4p we obtain 

r* ^ ,n P [*>*K+k < <] = 1 

(t,x)^(t ,K) 

and the result follows. □ 

Next, we deal with functions P(t,x) = F(t, e x ), and so it is natural to consider the domain Dq = 
(O,/0) x (fc, In if) for some positive k < In K (which is, essentially, the same domain after the space 
coordinate change) and let 

(4.11) c = P A ,t(P>K)= I™ PA,t(fi,x). 

x— »log K 

Let v(t,x) be a function solving the equation + Aw(i,x) = with A defined by Q3.6P and satisfying 
the boundary conditions 

(4.12) u(i,miT) = c , v(t,k) = P t (t,k) for < t < /3 and v(/3,x) = P t (/3,x) for k < x < hxK. 

Since these boundary conditions are continuous then (see [3]) they are satisfied by a unique solution in 
C 1,2 [D] of the above equation. Let w(t,x) be a function on Do such that 

(4.13) P t (t,x) = w(t,x)+v(t,x) V(t, x) G Do \ ((3, In K). 

Thus, w(t,x) G C 1,2 [-D'] and it satisfies the same parabolic equation in Do as ^r(t,x) and v(t,x). Its 
boundary values are 

(4.14) w(t,1nK) = -c , w(i,fc) = for < t < /3 and w(0,x) = for k < x < In if. 

From the continuity of v(t, x) on Do we see that it is bounded there and since ^ is also bounded there 
we obtain the same result for the function w as for v. Hence, 

(4.15) w(t, x),v(t, x) G C h2 [D ] n L°°[D ]. 

4.2. Integrability of wt(t,x) and w x (t,x). Now we will analyze the function w(t,x). Let z\ u ' x ^ — 
(u + t, Xf) be the diffusion process in the plane whose infinitesimal generator is equal to Li = + A 
on the space of C 2 functions. For each e > define D E — (0,/3 — e) x [k + e,\nK — s). Let T e be the 
parabolic boundary of D e . For every e > which is sufficiently small we can find a smooth function 
w(t, x) with compact support on the plane such that in D e it is equal to w(t, x). By the Dynkin formula 
we obtain that for every (u, x) G D e , 

MT e ) 

(4.16) E [U)J ;][i(j(^ T ( re ))] = w(u,x) +E[„ iX] [ / L 1 w(Z s )ds] 
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where t(Q) denotes the arrival time to Q by the process z\ u ' x K Note that since w(t, x) = w(t,x) for 
(t, x) G D e we can replace w by w in the above formula and since z\™' x ^ G -Do for s < r we obtain that 
Jjw(Y s ) — 0. It follows that for every e > 0, 

(4.17) w(u,x) = E [UiX] [w(Z T ^)]. 

Now fix (u, x) G Dq and a continuous path loq. Let £ = {zj^d^ • ) (wo)} no < n C Dq where uq is such that 
(u,x) G D i . The sequence of times {T(ri)(w)}„>„ is non decreasing with respect to n and so it has 
a limit p <T. Let 7 be an accumulation point in £, i.e. lim^oo Z^ x ^ )(^o) = 7 for some subsequence 

nk- Define d(y,To) = mf{\y — x\ : x G To} and note that this function is continuous on Do and it is 
if and only if y G To- Since d(Y^ x ^ ^(ujo)^o) < ^ for each k we conclude that 7 G To and since 

hnife^oo r(r^_) = p it follows that z\?' x \uq) = lim^oo Z^,' x \ Juj ) = 7. Hence, t(T ) (uj q ) = r. By the 

rtf, \ — ! 

definition w(t,x) is continuous except at the point (f3 7 lnK) but because P^z] [Z T (r ) = (/^ m -^0] = 
for every (u, x) G Dq we can ignore paths that reach the point (/3, In iT), and so 

(4.18) limHzgl^tzg,) P M a.s. 

4.4. Corollary. For every (t,x) G -Do, 

w{t,x) = E [ttX] [w(Z TiTo) )] = -cE [t , 2;] [I {r(rol)<r(ro2uro3)} ] = -cP [t)a .][r(r i) < r(r 02 U r 03 )] 

where T i = {(f.lnir) : < t < 0}, r 02 = {(fi,x) : k < x < \nK} and T a3 = {(t,k) : < t < /?}. 

Proof. From (|4. 1 5|) we know that the function w(t, x) is bounded and so we can use the Lebesgue bounded 
convergence theorem and from the boundary conditions on w(t, x) it follows that 

w(t,x) = lmLE [M] HZ r(re) )] = E [M] [lhn w(Z t{Tc) )\ = E [tM [w{Z T{To) )] 

which gives the first equality of the corollary while the second equality follows from (|4.14p and the third 
equality is obvious. □ 

Let (t,x), {f ,x) G -Do and assume that t < t' . Then it is not difficult to understand that 

p [t<x] [r(T 01 ) <r(r 02 ur 03 )] >p [t , fX] [r(r i) <r(r 02 ur 03 )], 

and so w(t, x) is nonincreasing in t for every x which implies that 

(4.19) —(t,x)>0 V(t,x) G D . 
It is also easy to see that for < t < T and < x < x' < 1, 

P [M] [r(r i) <T(r 02 uro 3 )] < p [Ml [r(r i) < r(r 02 ur 03 )], 

and so 

(4.20) V ('- I ) 6fl »' 

4.5. Lemma. The functions w t and w x are in L 1 [D ] . 

Proof. We will use (|4.19[) in order to prove the result for wt(t,x). The case of w x (t,x) can be proven 
similarly be using (|4.20j) . Using (|4.14j) . (|4.19j) and the continuity of w(0,x) on {0} x [fc,lniT] we obtain 
that 

Id„ \w\dtdx = f k nK f P ^dtdx = lim^o £ K ~' j* ^dtdx 
= lim e ^o f k ^ K £ (w(/3, x) — w(0, x))dx — — lim £ _>o Ji? S w (0, x)dx — — f. w(0, x)dx < oo. 
Using (I4.20P in place of (14. 19|) the proof of integrability of w x is similar. □ 
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4.3. Integrability of v t (t,x) and v x (t, x). We continue this section by analyzing the function v(t,x) 
solving the equation Liu = with the boundary conditions given by (|4.12l) . Let C 1:2 [Dq] be the set of 
all functions which have one derivative in t and two derivatives in x both uniformly continuous in Dq. 

4.6. Lemma. There exist a function z(t,x) G C 1:2 [l)o] such that 

(4.21) z(t,x) =v(t,x) v(f,aOer . 

Proof. Recall that P A j(T, x) = Pt(T, x) for k < x <\nK and note that the functions P Ajt (T, a;), P At t(t, x) 
and P t (t,k) as function of (t,x) belong to the space C 1,2 [Dq]. Set 

~z(t,x) = ^ ~ X k (P t (t, k) + P Att (T, x) - P A , t (T, k)) + ^x* k p A,t{t, x). 

Then z(t, x) G C 1,2 [Z)] since it is a linear combination of functions from this space. We also have 
^M) = l £%E%(Pt(t,k) + P A , t (T,k) - P Ajt (T,k)) =P t (t,x) V0 < t < (3 
z(t,\nK) = P Ait (t,]nK) when < t < (3 and for all k < x < InK, 
~ Z ^ T ^) = 1 iw^( p t(T,k)+P A>t (T,x) - P A>t (T,k)) + ^^P A ^ t (T,x) = P A , t (T,x). 
Thus, we obtain 

(4.22) z(t,x) = ^Z£i{t,x) + J£JL.z{T,z) € C^[D]. 
Since 

z(t,k) = z(t,k) = P t (t,x), z(t,hxK) = z(T,\nK) = P A}t (T,lnK) = c and z(T,x) = S{T,x) = P At (T,x) 
it follows that 

(4.23) z(t, x) = v(t, x) V(t, x) G T. 

□ 

Next, define f{t,x) = —~Lz(t,x). From Lemma 14.61 we obtain that f(x,t) is bounded in Dq and so it 
belongs to L p [Dq] for every 1 < p < oo. Set v(t, x) = v(t, x) — z(x, t) and observe that 

liv(t, x) = f{t, x) and v(t,x) = V(f, i) £ Tq. 

We conclude that the function v(t,x) is the unique solution of the following problem (see [T]). 

4.7. Theorem. Let 1 < p < oo then for any f(t,x) G L P [D ] there exists a unique function v such that 

(i) v G L»[0, T; W^iO, 1)] n W[0, T; W^ p (0, 1)}, 
(U) o*e L p [Dq]i 

(Hi) Lv{t,x) — f(t,x) for every (t,x) G Dq, 
(iv) v\ To = 0. 

From assertions (i) and (ii) of Theorem 14.71 we obtain that the functions v x (t, x) and Vt(t, x) are both 
in L P [D] for every < p < oo and since z(t,x) G C 1,2 [Dq] we obtain the following. 

4.8. Corollary. For every 1 < p < oo the functions Vt(t,x) and v x (t,x) belong to the space L p [Dq]. 
We can now summarize most of the results of this section as follows. 

4.9. Proposition. Let s(f3) < k < \nK < k' and define 

D Q = (0,(3) x (k,\nK) and D' = (0,(3) x (\nK,k'). 

Then the function Pt(t,x) is continuous at every point in the domain Dq \ {((3,hiK)}, and there exist 
two functions w(t,x) and v(t,x) on Dq such that 

(4.24) P t (t,x) = w(t,x) +v(t,x) for every (t, x) G Dq \ {((3, In K)}, 
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(4.25) w(t,x),v(t,x) G C 1 ' 2 (D )UL oo [D ] 

and both functions are solutions of the parabolic equation L^u — 0. Furthermore, w(t,x) is continuous in 
D and it satisfies 

(4.26) w(t,\nK) = P Att (J3,\nK) andw(t,b) = when < t < /3, 

w([3, x) = when k < x < In K 

and 

(4.27) w t (t,x),w x (t,x) € L^D]. 
Finally, v(t,x) G C(D) and for every 1 < p < oo, 

(4.28) 

XTie same decomposition of Pt(t,x) with the same properties holds true in the domain D f . 

Proof. Taking the same functions v and w as in (|4.13p we see that (I4.25[) is actually the same as (|4.15[) 
and the fact that both v and w are solution of Liu = is clear from their definitions. Next we see that 
(14.261) is the same as (|4.14p . that (|4.2T[) is the same as Lemma \4. 51 and that (|4. 28[) is, in fact, Corollary 
14.81 Observe that we did not use in this section the fact that k < In K so all the proofs are also applicable 



(4.30) 



to the case k' > \nK and the domain D' . □ 

From (f4~24|) . (|4T27|) and ([OS} we obtain the following. 
4.10. Corollary. Let D = {(t, x) : < t < (3, k-fit<x<\nK~ /it} and 
(4.29) u(t,x) = e^ rt P(t,x + fit). 

Then 

~dT 2 

4.4. Price function when initial stock price is large. Let F(t,x), P(t,x) and u(t,x) be as above. 
Recall that in the domain (0,T) x (InK, oo) the function P(t,x) satisfies the equation LiP = 0, it is 
continuous in the closure of [0,T] x [\nK, oo) and P(T,x) = (K — e x ) + = for x > In A'. Define 

(4.31) v(t,x) =u(T-t, -^=x + In K +\n\T) 

v2 

where u is given by (|4.29p and set G — (0,T) x (0,oo). It follows from Proposition 14.91 that 

(1) v(t,x) £ C U [G]UC[G], 

(2) v xx (t,x) = v t (t,x) for every (t,x) G G, 

(3) v(t, 0) = u{T - t, In K + 2\(j,\T) is continuous, 

(4) v(0, x) = for every < t < T, 

(5) v(0,x) is bounded (since P(t,x) is). 

Since a bounded solution of the heat equation in G is unique (see [3]) then for every (t, x) G G, 

f* dK 1 x " 

(4.32) v(t,x) = -2 — (t - T,x)v(r,0)dT where K(t,x) 



;uh: so 



9x V47rf 
xe 4 (*--> u(r, 0)dr 



47 7 (t-r)3/ 2 • 
Differentiating u we obtain polynomials Qk,n(s,x) such that for all k,n G N, 



d n td k x 



(t,x)= / g„, fc ((i-r)- 1 / 2 ,x)e3(^) U (T,0)dT 
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If N is large enough and c > then Qk,n((t — r) -1 / 2 , x) is a polynomial in (t — r) 1 / 2 and 1/x and 

it is bounded on (0,T) x (c, oo). Since sup y > a y N e~ v < oo for any JVeNwe can set y = 4 (f_ T ) deriving 
that for any N eN and (t, x) e (0, T) x (c, oo), 

< (ir/o ( i£ #On, fc (^(i-r)- 1 / 2 ))y Jv e-^(r,0)dT< ^ 
For some C = C(N) > 0. Hence, the following results hold true. 

4.11. Corollary. For any k, n positive integers k, n and c > 0, 

5^) S£ > [( 0,r)x(e,=o,]. 

4.12. Corollary. Let &(\nK + \n\T) < k' and G = {(t,x) : 0<t< (3, k' - (J, < x < oo}. Then 

^(t,x)EL^G}. 

5. Proof of main theorem 
We split the proof into two cases for x < In K and for x > In A". 

5.1. Case x < In K. We begin by proving the upper bound in (12.11[) . Since the option holder can exercise 
at time it is clear from the definition of P(t, x) in (|2.5[) that P(t, x) > i/j(x) f° r every x > 0. Furthermore, 
by Proposition 13 . 1 f iv) for each fixed t the function P(t, x) as a function of x is nonincreasing. Therefore, 
P(t,x) > P(t,\nK) = S when x < \nK. From the definition ([2.7J1 of the stopping time cr") it is not 
difficult to see that in the present case when < T, 

x + fia (n) + nB a{n) < In K, 

and so 

(5.1) P(^ n \x + f m ( ^ + K B ffW ) >S. 
Hence for every r £ 

T (n) 

we obtain, 

(5.2) E[e-" ACT< "' (^(z + /rf + «4 n) )I T < CT (») + SItMkt)] 

< E[e— a«tW ( P(r A CT (n) ; x + MT A (j( n ) + kB M w )] = E[u(t A ff W , lJ n A » w )]). 
By Proposition [ 



T'ItAit'" 1 ) 

(5.3) E[u(t A .("U^U)]) = «M + E[ E Vu ^ l ) h ' X $-i) h )\ 

3=1 

where, as before, x) = e~ rt P(t, x + fit). Taking the sup with respect to all r <E in the inequality 
(|5.2[) and using the fact that u(0, x) = P(0, x) we obtain that 



h" 1 (TAo-(" ) ) 

(5.4) pW^-PfO,!)^ sup E[ T Vu((j-l)h,X^)}. 

Thus, in order to bound Pj (x) — P(0, x) from the above it suffices to find an upper bound of the right 
hand in ([57 
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Next, we split the domain [0, T] xl into three parts 
(5.5) C = {(t,x) E [0,T - h] : fit + x > s(t + h) + \fi\h+Ky/h}, 

S = {(*, x) E [0,T-h] : fit + x< s(t) - \fi\h - kVK} and 
B = {(t,x) E [0,T- h] X M : s(t) - \fi\h- kVK< fit + x < s(t + h) + \fi\h + K\fh}. 



In order to estimate the right hand side of (|5.4[) we split it into three parts according to the domains C, 
S and B, i.e. 

(5.6) E[Ej:; (TACTW) 2M(j - l)h,X^_ l)h )\ = EEj:;<™ W) %((i - l)h, 

+EEj:; (TA<r(n)) ^((j - i)^^g , 2 1)fc )%-i) fcJ x w _ 1)fc)6 B]. 

By Proposition 13 . 1 T ii) after the time j3 the prices of the American and game put options coincide which 
enables us to conclude that u(t,x) = e~ rt PA{t,x + fit) for t > (3 and that the sets C t >p — {(t,x) E C : 
t > &t>p — {(*) x) £ S : t > (3} and fit>p = {(t, x) E B : t > (3} are the same as the corresponding- 
parts of the domains C, S and B introduced in [T7] for the case of American put options. Therefore, we 
can use the following results from Sections 4.2 and 4.3 in |17j . 

5.1. Proposition. There exists a constant C > such that for every t E 

T (n) 

(r/h)Vk p r 

(5.7) E[ £ ^h((i-l)^^i)JI%-i)^_ 1))l ) e d<C'(^) 4/5 ! 
where kp = min{/c : kh > /?}, and 

(5-8) E[ E ^((i-l)^^;2i)J%-i)^ ( ,_ 1)h)eB ]<^74- 



Observe also that P(t,x) — K — e x in the domain S, and so we can use there Lemma 2 from Section 
4 of [17]. 

5.2. Lemma. i*br every (t,x) E S we have T>u(t,x) < 0, and so 

i=i 

Thus, for an upper bound of the right side of (|5.4[) we can ignore the second term in the right hand 
side of (|5.6[) and estimate only two remaining terms starting with the first term in the right hand side of 



5.3. Proposition. There is a constant C > such that for all n E N, 

/i-'fTAir'"') 

(5-9) E[ £ |2M(7 - l^^x^^-i^-^Jec] < Cn- 3 / 4 . 

Proof. We have 

(5-10) EEf:; (TAfftB)) l^((j - l)h,X$± 1)h )\I (U -i )h ^ 1)k )ed 
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Proposition 15.11 provides a bound for the second term in the right hand side of (15.101) , and so it remains 
to deal only with the first term there. Note that if jh < o-W A fi n ) and (jh,xffi) G C then 

c x (j) = s(jh) - fijh + nVh < X$ <\aK- 2nVh - fijh = c 2 (j) 
where the equalities above are just definitions of ci and &2- Observe also that since x < \nK and jh < cn n ) 

(r 



then by the definition of the stopping times a^> the process + fijh does not exceed In K — 2Ky/h. 



By Proposition 

I f K V Q2 U 

(5.11) T>u(t,x) = - dy dz(z——(t + y 2 ,x + z) + S(u)(t + y 2 ,x + z)). 

kJ J_ Ky dtdx 

Relying on the same computation as in Section 4 of [17] we see that for (t, x) € C and x < InK — \fi\h ■ 

(5.12) \Vu{t lX )\<^— ds dz\^(s,z)\. 

Thus, for < j < kp, 

^(|^(^,^)|I {07iiX <„ ))eC}n{j7i<(T( „ )} ) < \Du(jh,y)\dP x Jy) 

< ^^ +hds ^ h ms,z)\dz)dV Xjh {y) 

From (|3 .4|) we see that there is a constant C > independent of j and n such that 

P[\X^-z\< K Vh]<^S=. 

Hence, for jh < a^ n \ 



n\Vu(jh,X^ jhx ^ < %f<r >n ds^^dz\^z)\^= 

ds 



2k Jjh Jh(i+X) J ci(])-KVh 1 9t 2 v ' n 



— n Jj/i Vs Jc 1 (j)-KVh 1 <9* v ' yl 

Define 

Ci(t) = s(t) — fit, 02(f) = In K — fit — nVh 

where s(t) = ln(&(t)) is the free boundary of the option holder and b(t) was introduced at the beginning 
of Section |3l Observe that for every j and any jh < s < (j + l)h, 

ci(i) - nVh > ci(s), c 2 (j) + Kv/i < c 2 (s). 
Summing up the above estimates we obtain 

= t + Ht(lf h f s Ltl z)\ + f?- h % C^dz\^(s, z)\) 
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where the term ^ comes from the first term E|2?u(0, x)\ of the sum which can be estimated easily using 
the fact that Ut(t,x) and u xx (t,x) are bounded for small t. 

Let G = {(t,x) : < t < (3 : ci(t) < x < \nK - fit} and note that G C E U E a where E and E a are 
defined in Corollaries 14. 101 and 13.141 which imply that §jt(s, z) g L l [F}. Hence, 

(5.14) / ^= d z\^( S ,z)\<C in ^ ds dzl^^z^KCn 1 / 4 . 

Next, we estimate the first integral in brackets in the right hand side of (15 . 13[) . Let s(/3) < k < In if, k' = 
XnK ^ k and split the integral in question as follows 

(5-15) lf h ^CSdz\^(s,z)\ 

= if fs & *)l + If * I^ts *)\- 

From Corollary 13. 141 we know that the function ^r(s, z) is in L 2 [E], where 

E = {(s, z) : < s < T, a(t) < z < k' - fit} C E a 

(for an appropriate b < In K in the definition of E a ). Therefore we can use the Cauchy-Schwarz inequality 
to obtain 

(5-16) lf%Scjrdz\^{s,z)\ 



< 



(lf h ^L%» S dz)(ff h l^ S l&MM < Gin, 



Now we are left with the second integral in the right hand side of (15 . 15|) . We will show that there is a 
constant C > such that, 



(5.17) X n =\ BL\ dz\^(s,z)\<Cn^. 



/s 

-rt 



Recall that u(t, x) = e P(t, x + fit), and so 

: (i, x) = e~ rt (r 2 P(t, x + fit) - 2rP t {t, x + fit) - 2rfiP x (t, x + fit)) 



d 2 u, 



dt 2 

+e- rt (fi 2 P xx (t, x + fit) + 2fiP xt {t, x + fit) + P tt {t, x + fit)). 

Observe that the functions P(t,x), P x (t,x), Pt(t,x) and P xx (t, x) are all bounded for small t. Indeed, 
P < K + 5 while Pt is bounded in the domain of integration in (|5.17|) for small h in view of (|4.13[) , (|4.15[) , 
(14.241) and (|4.25p . Next, P x is bounded by Theorem 8.1 from [TJ. Finally, P xx is bounded since in the 
domain in question P and its first derivatives are bounded and P satisfies the equation ( + A)P = 
(see (|3.8[1 ). Therefore, we can write 

r^h fji r c 2(t) 

(5.18) l n < — dx(\2fie- rt P tx {t,x + iJ,t)\ + \e- rt Pu(t,x + iJ,t)\)+C l , 

Jh vtJk'-iit 

for some constant G\ > independent of n. Recall that for (x,t) E D = (0,/3) x (fc, In if) by Proposition 
14.91 Pt{t, x) = v(t, x) + w(t, x) where v t and v x belong to L 2 [D]. Hence, expressing P tx and P u via v t ,w t 
and v x ,w x we can estimate the integral (I5.18[) containing vt and v x by means of the Cauchy-Schwarz 
inequality as it was done in (15.16[) . Replacing these integrals by C2 Inn we obtain 



rVh fa rc 2 (t) 

Zn< I —j=\ dx(\2fie~ rt w x (t,x + fit)\ + \e~ rt w t (t,x + fit)\) +C 2 lnn + Ci. 
Vt Jk'-fit 
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By (|4.19l) and (I4.20[) the functions w t (t,x) and w x (t,x) do not change signs in D, and so it follows that 
(5.19) * J^ t (\2fie- rt w x (t, x + fit)\ + \e- rt w t (t, x + fit)\)dx 

ff h fj^l t ^w t (t,x + nt)dx. 



jf h ^2\^f: K -^ Wx (t,x)dx 



By Proposition 14.91 w(x,t) is bounded on D, and so the contribution of the first integral in the right 
hand side of (|5 . 19[) is bounded by a constant and it remains to estimate only the second integral there. 
Next, we will need a more explicit representation of the function w. Let 



(5.20) 



z{t, x) 



t w(t 1 x + {.it). 



Then in the domain E — {(t, x), < t < f3, k — fit < x < In K — fit}, 

k 2 

— z xx (t,x) + z t (t,x) = 0. 

Define 
(5.21) 

and let 



z(t,x) = z(T-t,—x) 



E = {(t, x) : < t < T, -5— i ^ < x < — — }. 

K K 

In the domain E the function z(t,x) satisfies the heat equation 

z xx (t,x) = z t (t,x). 

If we let 



(5.22) 



^ )= V2( fc -M(T-*)) ; d2{t)= V2(lnK-,(T-t)) 



then from the boundary values of w(t, x) we obtain 

z(0,x) = for di(0) < x < d 2 (0), z(t, Sl (t)) = and z(t, s 2 (t)) = e" r(T_ * ) for < t < T. 

Note that z(t,x) is a bounded continuous function on the boundaries (t,di(t)), i = 1,2 , < t < T of 
E. Hence, by Chapter 14 of [3] we can represent z(t,x) in the form 



(5.23) z(t,x) 
where H(t, x) - 



* dH 
o dx 



(x-d 1 {T),t-T)(j) 1 (T)dT + 



dH 

dx 



(x - d 2 (r),t - T)4> 2 (T)dT 



4nt 



is the fundamental solution and the functions 4>i(t), i — 1,2 are bounded 



continuous on the interval (0,T]. From the definition of z we see that 

Zt(t, x) — —re~ rt w(t, x + fit) + e~ rt wt(t, x + fit). 
Since w(t, x) is bounded then for some constant C\ > independent of n, 



rVh dt ,c 2 (t) . vfl ^ p W 

Jh y/tJk'-nt Jh vtJk'-iit 



Vh dt ,c 2 (t) 



z t (t, x)dx 



From the representation (|5.23p of z(t, x) we obtain that 
(5.24) 



< 



V2 
72 



^ S S /o T ~* H (* - * W, T - t - r)4>i (r)drdx 

ft lip-,*) i So~* %{*-Mr),T-t- r)0 2 (r)drdx 
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Observe that as long as we keep x or t away from the function H (x, t) is smooth and it has bounded 
derivatives with bounds depending on the range of t, x and their distance from zero. Next, if x satisfies 

fey fey fc\ 

— ik' -/Mt)<x< —c 2 (t) = — QnK -fit - nVh) 

K K K 



then 

/2 

k'-k < — (k' -k + u(T -t-r)) < x-diM for < r < T - t. 

Since k' > k we see that x — d\(r) stays away from on the entire interval (0, T — t\. It follows from the 
above that the function 

rT—t qjj 

$l(t,x)= / {x-dl{T),T-t-T)<h.(T)dT 



has bounded derivatives with respect to t with bounds independent of n in the region {(t, x) : h < t < 
y/h, —(k' — fit) < x^-C2(t)}. We conclude that the hrst integral in the right hand side of (I5.24j) is 
bounded from above by a constant independent of n and it remains to estimate the second integral there. 
Set 

r T-t 



*2(t,a;) 



dH 

dx 



(x - d 2 {r),T - t — T)4> 2 {T)dT. 



We see that if 



then 



V2 \/2 
x < — c 2 (t) = — QnK - fit - Vh), 



x - daM = x-—(\nK - fi(T - t)) < —(u(T - t - t) - Vh). 

K K 

In this case x — cfafr) can be zero when t£ [0,T — t] but this can happen only for a r that which is at 
least /i~ 1 Vh apart from T — t. Thus, the function $2 is smooth with a bounded uniformly continuous 
derivative with respect to t though this bound may depend on n. Nevertheless, we still have the following 



V2 

— K 

~ V2 



ft IS £ ( Jo'' ^-d 2 (T),T-t-T)<f> 2 (r)dr) dx 



V2 



-t-T 



< 



V2 



-H(^(k' - fit)-d 2 (r),T-t- T))0 2 (r)rfr) 



+— 

+ \/2 



Sf fti{^ H ^ k ' -\nK + fi(T-t-r),T-t- r)Mr)dr 



We see that in the second term in the right hand side k' — In K + fi(T — t — t) can take on the value for 
t £ (0, T — t] but then r is at least c = 1 1 A;' — In K | apart from T — t and now the separation constant 
c does not depend on n. Thus, we can bound the second term there from above by a constant and it 
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remains to estimate the first term which we do as follows 

ffftiir* H{^c 2 (t) - d 2 (r),T - t - r)Mr)dr 



I — — 



n rVh dt fT-t 
^ C 2 J h Tt J 

n p</h dt i-T-t 
+ C 2 J h — t Jo 



(T-t-r) 3 / 2 eX P ' 



4(T-*-t) 



-)Mr) 



f2 r 



■ r* dt r-T-t 



(T-t-r)!/2 ex P V 4(T-t-r) 

«xp ( - aprfey) exp - ^(T - t - r))Mr) 



eh 



(T-t-r)5/2 V 2{T-t-r)) V ^2 

where C2 > is a constant independent of n. Analyzing the integral with respect to r in the second term 
in the right hand side above by considering different possible values of T — t — r we conclude that this 
integral is bounded by a constant independent of n. Next we observe that | exp jpz(T— t— t))</>2(t)| 

is also bounded by a constant independent of n too. Hence, we obtain 

"Vh dt rT-t 



I<C 3 + C 3 J7 > l fj*- 1 ^P ( - ) dr 

+ °3 if" fl Jq 1 (T-t-r)^/ 2 eX P ( ~ 2(T-t-r) ) dT 



for a constant C3 > independent of n. Set p 



„,„ h . — r and note that ^f- 

2\1 —t — T) ar 



Vh 

" 4(T-i—r) 3 / 2 



and 



(It 



' 4(T-t-r) 



-pp. We proceed by changing variables arriving at 



i < Ci + Ci c ft ^ + ^ r * vV^V 



< C 4 + C 5 ^ ^ < d + C 5 2(l + 1 rr 7I )< Cqti 1 ^ 

for some constants C4, C5, Cg > independent of n and (|5.17p follows. Combining (|5.17[) and (|5 . 1 6[) we 
obtain from (|5.15[) that 

Jh Vs Jc 2 (s) dt 

a 



(5.25) 



Finally, Proposition IBTTl follows from (j5~25l) . ([5~T3)) and (J5TT4J). 

Next, we turn our attention to the domain B. First, we will prove the following result. 
5.4. Lemma. There exists a constant C > such that for all n G N. 

(5.26) E[ Vu ^ ~ ^^S-i^kU-^x^^B] < Cn-*'\ 

3 = 1 

Proof. Let B t< o( n ) and B t>(3 („) be the set of all points (t, x) € B such that t < /3^ n ' and t > f3^ n \ 
respectively. We split (|5.26l) according to these two regions, namely, 



E[E£ (TAff( " ,) ^((j-i)^^ 1) Ji(( J -i, V(J _ 1)JeB ] 



+ e [Ej=i (tA<T< " >) Vu{[j - l)h, Xi;U )h )ku-i)H,x (j _ m MB^ (n) 



-u 

r (n) 



Vu((j - l)h 1 X (j _ 1)h )l a ^ 1)htXu _ 1)h)le - Bt 



-(») 



By Proposition 15. II we have that for a constant C > independent of n, 



3 = 1 
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Thus, it remains to estimate only the first term in the right hand side. Let E — {(t,x) : < t < 
p( n ',a — [it < x <b — fit} where a < s(0) and s((3^ + h) + \fi\h + 2n^fh < b < In A'. For n large enough 
we can find such a b because s(t) is continuous and s(P^) < In K. We know from Corollary 13. f 41 that 
u(t,x) € H 2 [E]. Since C 2 [E] is dense in this space we can approximate u(t,x) by C 2 functions to get 
equality (|3.3j) of Proposition 13.31 for u(t,x), as well. Since Ut(t,x) + ^Uxxit, x) < in the domain E we 
obtain 

Vu{t, x) < i Sf h dy JZ V dz(zit(t + y*,x + z)) 
<lfy<iyI%dz\^(t + y\ X + Z )\ 



It follows that 



Y rt+h i-s(t+h)+\Vh-nt q2 u 



where A = \n\ + K. Hence, 



EEj:; (TAaW) ^(^i)ft,^ 1) , 1 )i{(0'-i)^- 1) ,))eB i< , W }] 

Here kp = > an d the term — is the contribution of T>u(0, Xq) = Vu{0 1 x) < — which holds true from 
by the definition of the operator T> and boundedness of Ut and u xx for small t. From Corollary 13.141 we 
see that there exists a constant G\ > such that 

(5.27) / \-^-(t,z)\ 2 dz < d when < t < , 

Ja Otdx 

This together with p.4[) . the Cauchy-Schwarz inequality and the inequality -h= > 1., , which is satisfied 
when j > 1 and jTi < r < 2j7i, yields that 

EE £W»>) _ JC ( W 1)h )I {(0 ._ 1)h , Xo _ 1)h))eBt<;3( „ )} ] 

From Proposition l3.4l and Lipschitz continuity of the function P(t, x)'mt< p( n > uniformly in x < In K 
(see Theorem 8.1 in [14]) we obtain that for some constant C3 > 0, 



|«(*i) ~ «(*2)| < y/\h -h\C 3 whenever < h,t 2 < P {n) ■ 

Hence, 

fc -1 (TA/3 (B) ) 



Ef 



c 4 



n 3/4 

for some constant C4 > independent of n. □ 



By combining the results of Lemma l5.2l . Proposition l5 . 31 and Lemma 15^41 together with (|5.6[) we obtain 
that the upper bound P^ (x) — P(0, a;) < ^374 for some constant C > independent of n and of x < In K. 
Next, we will obtain a lower bound for the approximation error P-J (%) — P(0, x) when x <liiK. Set 

(5.28) T {n) = inf{t : n[t/h]h + x[ n) < s([t/h]h + h) + \fi\h + nVh}. 
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By Proposition 

(5.29) E[«(rW A .W,^^)] = u(0,x) + EE^^ \Vu((j - l)h, 

Set a = a n = T — ^3 and let be denned by (|5.28p with s there replaced by the free boundary sa 
for the American put option (see Section 2.2 in [T7]). Define also — T ( - n H^ T ( n)+h<a -j + TIi T {n)+h> a \ 
and Til, — T^A^t (>») , u _ 1 + ri f (») , L ^ t We will rely on the following estimate from Section 4.5 in 

A-i®- -rl \T A +h<a\ {t a +h>a\ J ° 

Hz]. 

5.5. Lemma. There exists a constant C > independent of nN swc/i i/iat 
(5-30) lEMrtUJ^-e^^+^w)]! < -£3 

where UA(t,x) — e~ tr PA(t,x + fit) with Pa given by &2.3\) . 



5.6. Remark. Note that s A (t) = s(t) for (3 < t < T, and so t^ V (3 = V /3 and V /3 = T a n) V /?. 

From now on we assume that n is large enough so that < a. From the definition of P[ n \x) we 
have 

(5.31) P[ n \ X ) > E[e-^^ n) ty^T™ + ^))I {r (n)< a (n, } + ^(-) <T W } )]. 

Hence, if we prove that for some constant C > independent of n, 
(5.32) 

J = |E[u(r(") A ff ("US, Aff «„))] - Ele-^'^'"' (^(rf + xW)I {T w< ffW} + fl W<T „)]| < ' 



n 



then by (|5.3ip and (|5.29[) we could conclude that 

(5.33) - -^=<P\ n) {x)- P(Q,x). 
We split the left hand side of (|5.32j) into three parts 

(5.34) J = E[{u(tW A[3( n \xl% w ) -e- rrln) ^ {n) (iP(fiT^ A /3< n > 
+^)A^)))}V<"»<^»A^)}]+E[{ M (a(™),X^ ) )"e— in) 6}l {aM<rkn)} ] 

+E[{u(r<»UgO - e^T^ 

This equality is true since = Ta = A /3 on the set Tq™' < <j( n ) A {3^ < a. We begin with the last 
term. First note that on the set < T a n ^ < we have, in particular, < cr 1 ™) and so = T 
by Remark 15.61 In the case Tq™' > we have r« = and r^") = r^ 1 ' and so from Lemma \5. 51 we 
derive that 

|E[{«(rW , Xg) - e""* (^"' + )) }l {/3( „ )<e »<^) } ]l 
< |E[Hri"\xW)-e-^i(^ T W ))}]| < ^. 

Next, we deal with the first term in the right hand side of (f5T34|) where T a n) = r (n) < er (n) A /3 (n) . This 
means that before time (3^ the process X^ n ' is stopped near the boundary s(t) and 



[IT 



W + < s(^ (n) + ft) + lM|ft + <n/fc. 
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By the definition, u(tW,I^) = e~ rT(n) P{t^ , l'") + jurO 1 )). Thus, we have 

E[{«(rW A ^\X%) - e — w a/*» A + X^ w ))}I {t ^ W} ] 

= E[{e-- ( " , (P(r(«),X(»)+ M r("))-^(MT ( " ) +xW))}l {r w 
If jUT^ < s(tW) then P(t("',I (b) + M^ (n) ) - = so we can assume that 

(5.35) s(r^) < fiT^ + Xfy < s(r (n) + h) + \fi\h + aVh. 
To continue we need the following lemma. 

5.7. Lemma. There is a constant C > independent of n such that for every point (t,x) satisfying 
s(t) < fit + x < s(t + h) + \fi\h + aVh andO<t < f3 (n \ 

C 

\P(t, /it + x) - ip(nt + x)\ < -. 

n 

Proof. The function P(t, x) is Lipschitz continuous when t < (3 and < x < fi(3 + \nK (see [2]), and so 

C 

\P(t,fj,t + x) - P(t + h, fit + x)\ < — 

n 

for some C > independent of n. If fit + x < s(t + h) then P(t + h, fit + x) — ip(fit + x) and we are done. 
Now assume that s(t + h) < fit + x < s(t + h) + where A = ^/n{\fi\h + <j\fh). 

From Corollary 13. 131 it follows that for every t <T and a < InK the function P xx (t,x) is continuous 
in x on the closed interval [s(i), a], so we can write 

A A 2 

Pit + h, fit + x) = P x {t + h, s(t + h))—= + P xx (t + h, s(t + h)) ha 

y?i In 

where a = a(h) satisfies lim^o(r) = 0. From the property of smooth fit (see [2]) it follows that 
P x (t + h, s(t + h)) — ip x (s(t + h)), and so for some C > independent of n, 

C 

\P(t + h, fit + x) - ip(ut + x)\< —. 

n 

□ 

Using (|5.35[) and the above lemma we obtain 

(5.36) |E[{u(r(") A ^"US^w) " e^'"^'"' (^(/jt'") A 

+^r<")A / 3(")))} I[ {ri" ) < -(™)A/3(™)}]l - n' 

Hence, we are done with the first term in the right hand side of (|5.34l) and it remains to estimate the 
second one. Since < < T the process X^ is stopped near the writer's boundary. Namely, we 
have 

h\K- \fi\h - aVh < fm {n) + X^ < InK. 

Since P(t,\nK) = S when t < (3, (3^ — f3 < h and P is Lipschitz continuous (see Theorem 8.1 of [Ti] ) 
we obtain that 



for some C > independent of n. Hence, 

(5-37) E[( ll (a(»),lg ) )-e- W i)l {aW<r , )} ] < ^ 

It follows that there exists C > independent of n such that for every x < InK, 
(5.38) <P$ n \x)-P(0,x). 
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Next, we will derive a lower bound for the second approximation function P^ n \x) denned by (|2.10p . 
still assuming that x < In K. According to (|5 . 29[) in order to obtain 

(5.39) f^-PtO,^-^. 
it suffices to show that 

(5-40) E[«(t<»> A oM,X$^)] - Pi n) {x) < 
We have 

(5-41) EKrWA.W.lW^Jj.pW^) < 

+ (^ n) + + *)V) <T W})] = E[{«(rW A 0<">,*$, A/jCn) ) - e— (n) ^ (n) (^r^ A 

+xS a ^)}I {t c^^ 

+ E[{u(rW,^ ) )-e-^ ) (^ M rW+x| ) ))}l ww<r , %ff( „^ 

Indeed, the first inequality is true since P^ (x) is defined as the sup on r € T (n) and we chose a specific 
one, i.e. t4 ■ The equality is true due to the same reason that 1(5.340 holds true. We see that the first 
term in the right hand side of (|5.41[) is the same as the first term in (|5.34l) and by (|5.36l) it is less then — 
for some constant C. The second term is nonpositive because for every (t, x) we have P(t, x) < tp( x ) + $ 
and u(t, x) = e~ rt P(t, /it + x) so we can just remove it from the equation. The last term is the same as 
the last term of (|5.34[) and from Lemma (|5.5I) we obtain that this term is less or equal than ^3 for an 
appropriate C. These arguments yield (15.40)) and hence (|5.39p . as well. For the upper bound we already 
know that p[ n \x) — P(0,x) < -^jj and from the definition of p[ n ^ and P^ it is not hard to see that 
|PjW — pj^\ < -£j=. It follows from above that there exist C > such that for every x < hiK, 

(5.42) -|<PW(,)-P(0,,)<-^. 

5.2. Case x >\rtK. We begin with the upper bound on P\ . We will show first that 

(5.43) P< n \x) - P(0,x) < sup E[ ]T Vuiij-l^xlf )]. 

The proof is similar to the proof of (|5.4p , we just have to show that for every r 6 T^ n ) , 

(5.44) P(rAaW )M rAaW+lW (n) ) 

> + X("))I {r < CT(n)} + (S - Ke(\f,\h + 2KVh))l {aW<T} . 

On the set r < cr'™) this inequality is clear since P(t,x) > ip(x). For the case u^ n ' < r observe that 
because x > In if we must have 

In K < na {n) + X^l < In K + \fi\h + 2kVK. 

By Theorem 8.1 in [TJ] the right derivative F x (t,K+) at K satisfies > F x (t,K+) > —1 for any t, and 
so < F(t,K) - F(t,K+ CX) < CX for each C > provided < A < A(C) is small enough. Assume 
< A < 1, then e A — 1 < Ae A < Ae. Hence, taking C — Ke we have 

P(t, In K) - P(t,\nK + A) = F(t, K) - F(t, Ke x ) < F(t, K)-F(t,K + KeX) < KeX. 
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Put A = \n\h + Ky/h then for crW < r and sufficiently large n, 

P(cr {n) , /itrW + Xfy ) > P(cr (n) , In if + A) > 5 - KeX. 

Hence, we obtain (|5.44l) which yields also (|5.4f I) . To bound the right hand side of (|5.43[) we split it 
similarly to the case x < In if (see (J5TSI) ) according to the three different regions C, B and S. Since our 
process starts at x > In if, if (Jj — l)h, € B for some j then this must happen after the time ft, 

and so we can use (|5.8[) . The part that belongs to the region S is non positive so we can ignore it, and 
so we will be left only with the region C. 

5.8. Lemma. For the discrete process x[ n ^ such that X^ — x > In K we have 
r'(TAiT w ) 

E[ Y, V \ u ^ - 1 )^^i)JI I {(0--i)'^o- 1 , ) J^}] < Cn-y\ 

3=1 

Proof. It suffices to show that 

fc/3A(CT < "'/'») 

(5.45) E[ Y 2? l«((j-l)^^i) h )|I{(( 3 --i)/.,x - 1)h )6C}] <Cn" 3 / 4 

3=1 

for some C > independent of n since after time f}( n ' we come back to the American option case. 
This is done in the same way as in Proposition 15.31 and so we provide only a sketch of the proof. Let 
c(s) = In if — fj,s + K-\/h then similarly to the proof of Proposition 15 . 31 we obtain 



Let ^(lnif + \n\T) < k' and split the integral in (|5.46l) into two parts 
(5-47) J,f$JS)«k|g?(«.*)l = 

Let £7 = {(s, z) : < s < /3, k! — [is < z < oo} then by Corollary 14. 121 we see that |jt(s,z) € L 2 [£], and 
so we obtain similarly to (I5.16|) that for some constant C > 0, 

p ds r°° , ,d 2 u, 



(5.48) / -= \ dz\ — (s,z)\<C\nn. 

Jh V s Jk'-ns Ct 

In the first integral in the right hand side of (|5.47p we do the same procedure as in (|5.13[) - (|5.1 Tl) relying 
on Proposition 14.91 and deriving that for some constant C > 0, 

(5.49) / ^ / dz\^(s,z)\<Cn^. 

Jh V* Jc(s) Ot Z 

Combining (|5.46p ~ (|5.49[) we obtain (|5.45p and complete the proof of the lemma. □ 

An estimate for the lower bound of p[ n \x) — P(0,x) when x > In if is done similarly to the case 
x < In if. As in that case we use the stopping time T^ n ' from (|5.28l) and from the above we see that 
(|5.29[) is true also for the case under consideration. We consider again defined before Lemma [5.51 
and similarly to (|5.30|) obtain that 

(5.50) |E[«(rir n) A ff W,X T (n) A|r(B) )] 

_ E[e -rrC»W(»> ^ (mt W + X^)I {rkn) ^ n)} + (S- Ke{MVh + vh))l {aln)<r p } ))]\. 
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In order to estimate (|5.50l) for i>lnJf we only need to split it into two parts, one for Ta < c'™- 1 and the 
other one for < ri . This it true in view of the fact that if we begin with x > ItlK and r^f 1 < 
then we must have (3^ < r« , and so we are back to the American option case and can use Lemma 15.51 
for this case. If cj(") < rt ] then the process is stopped near the seller's boundary and similarly to 
(|5.37p we can use the Lipschitz property of P to obtain, 

From here we can proceed similarly to the case of x < In K and obtain the lower bound for P^ proving 
(£10) for p[ n) . □ 

(n) 

Next, we turn to the second approximation function P 2 i still in the case of x > hiK. For the upper 



bound we use Lemma 1 5. 71 as in the case x < In K and proceed similarly to the proof of the upper bound 
for the first approximation function P^ . The proof of the lower bound is similar to the case x < In K 
and we obtain the result observing that if x > In K then P(t, x) < ip(x) + S = 5 for any t € [0,T]. □ 

6. Computations 

In this section we exhibit computations of price functions and free boundaries of game and American 
put options. All graphs of functions related to game put options were plotted using the approxima- 
tion function p^ 2000 - 1 (see (|2.10p V The graphs for the American put options were computed using the 
approximation function p^ 2000 ^ from [17] . 

Figure Q] shows both free boundaries of the holder and of the writer of a game put option and also 
the free boundary of the holder of an American put option corresponding to the option parameters 
K = 20, r = 0.02, k = 0.15, T — 0.5, S — 0.15. Here K is the strike of the option, r is the interest rate, 
k is the volatility, T is the time to maturity and S is the writer's cancelation penalty in the case of game 
option. 

In Figure [5] we plot the graphs of an American put option price function and of a game put option 
price functions with 5 = 1.0 and S = 1.5 while other parameters are K — 20, r = 0.02, k = 0.15, T = 10. 

Figure [3] shows the holder's free boundary for American and game put options where we use the same 
parameters as in |T]) adding also plots of free boundaries for the game put options with penalty values 
S = 0.3 and S = 0.5. 
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0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 

Figure 1. Free boundaries of American and game put options. 




Figure 2. The price functions of American and game put options. 
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Figure 3. Holder free boundaries of American and game put options. 
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